Understanding entangled spins in QED 
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The stability of two entangled spins dressed by electrons is studied by calculating the scattering 
phase shifts. The interaction between electrons is interpreted by fully relativistic QED and the 
screening effect is described phenomenologically in the Debye exponential form e~ ar . Our results 
show that if the (Einstein-Podolsky-Rosen-) EPR-type states are kept stable under the interaction 
of QED, the spatial wave function must be parity-dependent. The spin-singlet state s = and the 
polarized state -^=(| + — )— | — +)) along the z-axis give rise to two different kinds of phase shifts. 
Interestingly, the interaction between electrons in the spin-singlet pair is found to be attractive. 
Such an attraction could be very useful when we extract the entangled spins from superconductors. 
A mechanism to filter the entangled spins is also discussed. 
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Quantum entanglement has wide applications in many information processing protocols, such as quantum teleporta- 
tion [1], quantum communication and quantum computation [2], and the development of these quantum technologies 
depends strongly on both the properties of the entangled states and our capability to produce them. The strongly 
correlated quantum states were first obtained for photons in 1950 [3], and soon after Bohm expressed the EPR state 
with entangled spins [4]. Following these approach, Bell put forward his famous inequality [5] in 1964, and it was 
firmly demonstrated years later in experiments involving only massless photons [6]. Recently, the entangled spins 
t-H . dressed by relatively isolated electrons in solid-state materials appear to many researchers as a promising candidate 
' for the carrier of the EPR states that are required by quantum information processing [7] . 

The mechanisms for generating the entanglement of electron spins fall roughly into two types: one concerns specific 
systems in solid states [8], such as generating entangled electron spins via quantum dots or magnetic impurities, 
Oh! extracting Cooper pairs out of superconductors and etc.. The other employs special techniques by combining en- 
• tanglement with the effects of identity [9, 10] , special relativity [11-14], or electron-electron interaction within the 
entangled pair [15, 16, 17, 18]. Among the second type of mechanisms, of special interests are the space-spin entangle- 
ment transfer [10-13] and generating entanglement by dynamical scattering [17, 18]. However, serving as information 
carriers, the entangled spins are required to be stable enough when used to store and transport information. The 
question comes up is what conditions should be satisfied by the spatial wave function if the spin-singlet pair is to be 
kept stable under the interaction between electrons in solids. 

On the other hand, a difficulty encountered in entanglement generation is to overcome decoherence carsed by the 
; I ■ separation of an entangled pair [19] and the interaction between constituent particles [20]. Therefore, a full under- 
standing of the electron-electron interaction in the entangled pair is helpful in the design and control of entanglement. 
One category of theoretical mechanisms (thought experiment) to produce electron entanglement mainly cope with the 
interaction between electrons through the process of scattering. Previous studies either employ only an oversimplified 
Hamiltonian [15, 16], or analyze the scattering processes using quantum statistics or some elaborately designed inputs 
[18] instead of calculating the amplitude of scattering. In this paper we will describe the intra-pair spin interaction 
in a complete form covering all terms related to the spin operators and calculate the amplitude of scattering. Hence, 
knowledge gained in this work is fundamental and will be instructive in constructing any realistic models of generating 
entanglement of fermions. 

In what follows we focus our investigation on the stability of the EPR state r^(|TI) — lit)) by calculating the 
scattering phase shifts. The interaction between entangled electrons is described by fully relativistic QED. The 
method of phase shift has been widely used in many areas such as particle physics and nuclear physics, but it can only 
be applied to specific systems when interactions fall off more rapidly than coulomb potential 1/r [21]. In nature, as 
we know, it is hard to find pure coulomb potential. Therefore, if we are to examine the interaction between electrons 
that are entangled together, it is advisable to include the screening effect. Here, motivated by the idea of applying 
Debye theory [22] to the electrolyte, fluid, and dilute ions [23], we introduce a simple factor e~ ar (r is the distance 
between the two electrons) to phenomenologically illustrate this effect. In this picture, the total potential actually 
falls off more rapidly than coulomb potential. In solid, while the electrons are moving, the crystal lattice made up of 
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ions with positive charges and spins around the electrons will be distorted, yielding a screening effect. 

To understand the generation and evolution of the degree of entangled electron spins, Liu and Chen studied the 
interaction between the entangled electrons [16]. The authors analyzed the entanglement of two identical electrons with 
an interaction (neither covariant nor screened) interpreted by the nonrelativistic QED, and showed that the entangled 
spin-triplet states can evolve into states bearing no spin-entanglement, whereas the spin-singlet state remains stable 
in the scattering process. However, as shown in the following, our results suggest that to make the spin-singlet state 
stable under the electromagnetic interaction, some limitations on the angular momentum must be added to the spatial 
wave function. 

Here, we take the spin-singlet state of two electrons as initial and final scattering states to calculate the phase 
shifts. To make a comparative study, two methods of directly making total spin s = (Method A) [unpolarizing 
case] and polarizing the orientation of the constituent spins (Method B) are employed. To avoid confusion, we use 
-i=(||J,)— ||t)) to denotes general s — eigen state appeared in Method A and ^(| H — )— | — h)) the polarized 
state along the z-axis appeared in Method B. In quantum mechanics, the above two expressions usually both denote 
a spin-singlet state. In QED, however, since we calculate the unpolarized scattering process by averaging initial spins 
and summing over final spins, the results will be dependent on the polarization. So here, we use these two notations 
to distinguish two sorts of process (Method A and B). We will see that although the total spin for both processes 
are s = 0, other features are drastically different. Regarding entanglement, the two states are similar; but for the 
derivation of scattering amplitude, they differ from each other. 

For the evaluation of the electron-electron scattering phase shift with the total spin s = in the initial and final 
states by means of Method A, we need the amplitude of the scattering process M which has the following form [24] 
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u(p) is the Dirac spinor defined as J E 2 ^ { gj p ) > E = Vp 2 + m2 - Here the amplitude form is in fact an 

V V E+m J 

operator form, with the wave functions for spins determined by scattering characteristics. If the injected states and 
scattered states are not polarized, we perform the calculations using Method A, otherwise we use Method B. The 
indistinguishability between the entangled electrons is automatically satisfied by the amplitude. Since this amplitude 
is covariant, it is allowable to choose a special reference to simplify the formalism but at the same time leave the 
final matrix elements unchanged. Here the center-of-mass(CoM) reference system is used. Then the electrons' initial 
momenta satisfy pi = P2 = p, and the final one has P3 = P4 = q. For the elastic scattering process, the relation 
I p | = | q I is hold for the momenta. Substituting the Dirac spinors of the CoM into Eq. (1) leads to the explicit forms 
of Mi and M 2 [25]: 
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If the 7-matrix 7^ in Eq. (1) changes to 70, the first part of the amplitude reduces to the one that only the interaction 
of point charges is included (without the interaction of magnetic moments); and if furthermore the velocities of the 
two electrons are very low, the interaction can be approximately described by the classical Coulomb form i [26]. In 
the case 7^ — ► 70 Eq. (2) and Eq. (3) reduce respectively to 
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The screening effect can be phenomenologically incorporated into the propagators in Eq. (1) by introducing a factor 
e~ ar , where a is the inverse of Debye screening length [22] and r the distance between the two electrons. The factor 
can be related to the momentum propagator by applying the Fourier transformation 

/e — e*' ^ 3 k=-i-,. (6) 
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Obviously, if a — > the propagator will reduce to the original form. Now we can include the screening effect by using 
the propagator l/(fc 2 — a 2 ) instead of 1/fc 2 in Eq. (1). The phenomenal factor e~ ar in the propagator is effective for 
both scalar and vector potential, suggesting that the magnetic moments are screened in a similar way as charges. The 
screening effect of magnetic moments has been confirmed by Wilson et al in studying the Kondo effect [27] . Although 
the Kondo model has been studied intensively, its entanglement structure is still unclear [28]. Up to date no studies 
of scattering for entangled screening electrons have been reported. 

Now the phase shifts can be calculated under Born approximation, as we did in a previous study [29] 

5f = ~EkMj\{k). (7) 

Here E is the total energy of the two-electron system, k =| p |=| q | is the magnitude of the relative momenta p and 
q, and A4j\(k), with total angular momentum J and orbital angular momentum I, is the transition amplitude given 

by 
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where C^Jfi m are the Clebsch-Gordan coefficients, Y\ m (k) the spherical harmonic functions and Mfi(p, q; m s , m' s ) 

the matrix elements obtained directly from Eq. (1) by considering the total spins in the initial and final states: 
■Mfi(p, q; m s , m' s ) — (sm s ' \ Mfi \ sm s ). For special cases that s = or I = 0, i.e. without spin-orbit coupling, 
Eq.(8) reduces to [30] 



M d k ) = ^Zl dxPi{x)M }i {p,^x), (9) 

where Pi(x) is the Legendre polynomial with x = cos 9, and 9 is the angle between p and q. 

Now, let us evaluate the scattering phase shifts by means of Method B. The calculation of polarized amplitudes is 

analogous to that of a previous study [24], whereby the spinor u(p) in Eq. (1) is changed to J E ^ n ( J^ P ] x\ to 

v V E+m J 

include the spin states \\ [ A = 1,2 denote the two spin cigen states ^ q ^ an d ^1^1" -Here we assume that the 

momenta of initial electrons lie along the z-axis, and without losing generality, the spin directions lie along or opposite 
to the z-axis. The matrix elements from all possible combinations of polarized incident and scattered electrons have 
been listed in Ref. [24]. By imposing a proper transformation with respect to the defined polar angles of spins on 
these matrix elements, any other amplitudes with arbitrarily defined polarized orientation of the initial and final 
electrons can be obtained [24]. Expressing the {s = 0, s z = 0} state with the z-componcnt of spin angular momentum 
as ^=(| H — }— | — h)), the scattering matrix element between two {s — 7 s z — 0} states can be formally interpreted 

as -^«+- | -(-+ \)Mfi^{\ +-)- | -+)) = (-+ | M fi | -+) - (+- | M fi | -+} , where the two terms of RHS 
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are among the aforementioned list [24] . As there is no ambiguity in defining the related total angular momentum for 
partial waves with vanishing total spins, it can be postulated that Eq. (9) still works in this case if we replace the 
magnitudes with the polarized ones. 

The signs of the phase shifts can be determined by realizing that while the interaction between the electrons is 
interpreted as the classical Coulomb form [Eq. (4)], the interacting force of the S- wave is repulsive, and thus the 
phase shifts are negative. The resultant phase shifts of S-, P-, D-, and f-wave from Method A are plotted in Fig. 
1. The S- and D-wave phase shifts obtained with Method B are shown in Fig. 2 [31]. The P- and .F-wave shifts 
obtained with Method B vanish for the reason we will discuss below. The phase shifts will not change their signs in a 
wide region of a provided that a is lower than the electron mass. Generally, a larger a corresponds to smaller phase 
shifts. The dependence is illustrated in Figures 1 and 2. 

Phase shifts obtained with Method A (Fig.l) and B (Fig. 2) show much different features. The phase shifts from 
two methods for any given partial wave possess different signs. The sharp difference is presumably due to the use of 
different entangled spins in production as the scattering initial and final states. Although the two entanglement states 
all satisfy m s = 0, those obtained with Method A seem valid for all directions and those obtained with Method B only 
for z direction. The most salient common feature of the two sets of phase shifts is that they are both parity-dependent. 
Fig. 1 shows that the forces of S- and .D-wave are attractive and those of the P- and F-wave are strongly repulsive. 
We recall that according to the definition of parity, (—1)', / is the quantum number in spherical harmonics Yi m (9, ip), 
the S- and D-wave are even and the P- and F-wave are odd. In Fig. 2 we see that the S- and F>-wave are repulsive, 
the P- and F-wave, however, are all vanishing. The understanding of the parity-dependence may follow the fact that 
the even-parity spatial function of the two electrons combined with their s = antisymmetric spin wave function 
make up of the totally perfect antisymmetric wave function required by identical electrons. Hence, in Fig. 1 the 
states with the P- and F-wave spacial functions are forbidden by the strong repulsive force. In Fig. 2, however, the 
forbidden states of P- and F-wave are automatically removed by the special polarization in which the spin direction 
and its relationship to the momenta (spatial wave function) are defined simultaneously. From non-relativistic QED it 
is impossible to obtain the parity-dependence for the spatial wave function. 

It's clear that the interaction in Method A is attractive, and that in Method B is repulsive. The attractive force 
deserves more attention. Its order of magnitudes can be evaluated directly from the resultant phase shifts. With 
the assumption that the phase shifts are approximately in proportion to the transmitted momentum k in the region 
under concern, Eq. (7) yields dSi/dk w -2M | V | *j) « -2MV/(197) 3 . For a = 1, substituting the electron 
mass M = 5 x 10 5 eV and dSi/dk « 10~ 8 gives V ~ — 10~ 8 cV. For a smaller a, (e.g., 0.001), in low k energy region, 
V w -10~ 4 eV, on just the same order of magnitudes as the force of Cooper pair in a superconductor. The results 
might be heuristic in the development of the spintronic devices. The contributions of each term in Eq. (2) and Eq. 
(3) to the attractive and/or repulsive forces can be numerically determined in a straightforward way. The attraction 
mainly comes from the contribution of ov pi ui- P2 + &i- P2 02- pi, and the repulsion from o\- 02, ov Pi 02- pi and 
<ti • P2 02 • P2 • Contributions from the remaining terms such as the purely coulomb term p- are essentially cancelled 
out by the subtraction of two terms in Eq.(l). 

The characteristic of phase shifts obtained with Method A and B will not change if we take two-photon processes 
[16] and radiation corrections into account. For the two-photon processes, a similar substraction of Eq.(l) holds too, 
so the leading contribution of the processes is canceled out. The remaining terms will not change the signs of the phase 
shifts because, multiplied by the square of coupling constant, they are negligible to the next-leading-order terms of 
the tree level contribution. So, there is no need to consider the two-photon processes here. The radiative corrections 
can be done by replacing the masses and charges in scattering amplitudes with the effective ones [32] . They will not 
change the sign of the calculated phase shifts for there is no electron's propagator in the amplitude. To put it another 
way, the screening effect can be viewed as a part of the renormalization effect [33]. 

In light of the above results, we propose a mechanism to produce entangled pairs of massive spin- 1/2 particles. 
Electron pairs are first injected into a semiconductor [34] with aforementioned screening environment. In strong 
magnetic field spin-orbit coupling can be ignored, and then a particular relative space- wave- function of pairs, say, 
S'-wave or D-wetve can be filtered out at the beginning by adding a strong background magnetic field. Then with 
a electric field voltage V the pairs can be led to a square potential well with a thin square potential barrier in the 
middle which divides the well into two parts (see Fig. 3), the well and barrier can be formed by normal semiconductor 
layers [35]. Suppose that only one pair is allowed in a well for simplicity. The potential step of the well and the height 
of the dividing barrier are adjustable just as in practice. If the interaction is repulsive, the two electrons tend to be 
separated by the barrier and the possibility for either of them to tunnel through the barrier is low; whereas if the 
interaction is attractive as in the aforementioned spin-singlet state (Method A), the electrons tend to stay together 
in one side and the barrier is easy for them to tunnel through [36] . Thus the electrical conductivity of the well is very 
low in the former case and very high in the latter. In this way the pairs that make the well very conducting can be 
filtered out as the maximally entangled spins (in Method A). 

In summary, we have extensively examined the properties of the interacting entangled electrons in a fully relativistic 
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formalism with two different methods. Although the techniques used in this work to deal with the scattering processes 
are standard and straightforward in QED, the idea to apply them to electrons in solid-state is very meaningful. Also, 
the calculation of tensor forces required in this application is quite sophisticated and not straightforward. The parity- 
dependence of the phase shifts yielded from both approaches suggests that if the entanglement of spins is kept stable 
under covariant interaction, the selection of the spacial wave function will not be arbitrary. Furthermore, we find 
that the spin-singlet pair ^(|TI) — lit)) an d the polarized state -^(| H — )— | — h)) correspond to two different types 
of phase shifts. In the former the electron-electron interaction is attractive and in the latter it is repulsive. In our 
thought, it is very important to demonstrate by deliberate calculations that there is an attraction between entangled 
electrons in screening environment. The attraction of like-charged colloids confined between walls has been observed 
experimentally [37] , and the screening electron's attraction has since been supposed to exist by some researchers [38] . 
However, experimental evidence and theoretical basis for the latter were still lacking before our work report here. 
The attraction in principle can be tested by experiments in solid, and might be helpful in developing some devices of 
spintronics. 
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I. FIGURE CAPTIONS 

Fig. 1: The calculated electron-electron scattering phase shifts for S, P, D, F- waves from Method A. The dashed, 
solid, and dotted lines correspond to a = 0.1, 1 and 10, respectively. 

Fig. 2: The calculated electron-electron scattering phase shifts for S, P, D, F-waves from Method B. The dashed, 
solid, and dotted lines correspond to a = 0.1, 1 and 10, respectively. 

Fig. 3: A square potential well with a thin square potential barrier in the middle, and a confined electron pair with 
S- or £>-wave interactions. In spin singlet, the electrons attract each other and the device containing such wells will 
display high electrical conductivity; whereas in polarized case, the electron pairs will contribute little to conductance. 
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